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ABSTRACT 

We present 107 new epochs of optical monitoring data for the four brightest images of the 
gravitational lens SDSS J1004+4112 observed between October 2006 and June 2007. Combining 
this data with the previously obtained light curves, we determine the time delays between images 
A, B and C. We confirm our previous measurement finding that A leads B by AtBA = 40.6 ± 
1.8 days, and find that image C leads image A by Atqa — 821.6±2.1 days. The lower limit on the 
remaining delay is that image D lags image A by Ar^c > 1250 days. Based on the microlensing 
of images A and B we estimate that the accretion disk size at a rest wavelength of 2300A is 
1.q14.8±o.3 cm for a disk inclination of cosz — 1/2, which is consistent with the microlensing disk 
size-black hole mass correlation function given our estimate of the black hole mass from the Mgll 
line width of log Mbh/Mq = 8.44 ±0.14. The long delays allow us to fill in the seasonal gaps and 
assemble a continuous, densely sampled light curve spanning 5.7 years whose variability implies 
a structure function with a logarithmic slope of 7 = 0.35 ± 0.02. As C is the leading image, 
sharp features in the C light curve can be intensively studied 2.3 years later in the A/B pair, 
potentially allowing detailed reverberation mapping studies of a quasar at minimal cost. 
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The quasar SDSS J1004+4112 at z s = 1.734 is 
split into five im ages by an intervening galaxy clus- 
ter a t z,. = 0.68 (llnada et al. 20031 : llnada et al. 20051 



Qguri et al. 20041) . With a maximum image 
separation of 14'.'62, it is a rare example of 
a quasar gravitat i onally lensed by a cluster 
(jWambsganss 20031: llnada et al. 20061) . One of 
the most interesting applications of this system is 
to use the time delays between the lensed im- 



ages to study the structure of the cluster. If 
we assume the Hubble constant is known, then 
the delays break the primary model degener- 
acy of lensing studies (the "mass sheet degener- 
acy"), and the delay ratios constrain the struc- 
ture even if the Hubble constant is unknown. 
After its discovery, several groups modeled the 
expected time delays in SDSS J1004+4112 and 
their de pendence on the mean mass profile of the 



tncir depend ence on tne mean mass pronie or tne 
cluster dKawano fc Qguri 20061: lOguri et al. 2004 : 
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Williams 20041) . When we measured the shortest 
delay in the system, between images A and B, 
we found a longer delay than predicted by the 
models (Fohlmeister et al. 2007, hereafter Paper 
I) where the discrepancy probably arose because 
the models included the cD galaxy and the cluster 
halo but neglected the significant perturbations 
from the member galaxies. As we measure the 
longer delays, where the cluster potential should 
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be relatively more important than for the merg- 
ing A/B image pair, we would not expect cluster 
substructures to play as important a role. 

We also expect this lens to have a fairly short 
time scale for microlensing variability created by 
stars either in the intracluster medium or in galax- 
ies near the images. The internal velocities of a 
cluster are much higher than in a galaxy (700 km/s 
versus 200 km/s), and SDSS J1004+4112's posi- 
tion on the sky is almost orthogonal to the CMB 
dipole (Kogut et al. 1993), giving the observer 
a projected motion on the lens plane of almost 
300 km/s. In Paper I, we detected microlensing 
of the continuum emission of the A/B images in 
Paper I and there is also e vidence for microlens- 
ing of the CIV br oa d line (iRichards et al. 2004 ; 
Lamer et al. 20061 Gomez- Alvarez et al. 20o"^K 



Once we have measured the time delays we can 
remove the intrinsic quasar variability and use the 
microlensing variability to estimate the mean stel- 
lar mass and stellar surface density, the transverse 
velocities, and the struc tu re of the quasar source 
dGil-Merino et al. 20051: iMortonson et al. 20051 
IPoindexter et al. 20071 : Morgan et al. 2007l )~ 

Finally, we note that SDSS J1004+4112 could 
be an ideal laboratory for studying correlations 
in the intrinsic variability of quasars. With, im- 
age C leading images A and B by 2.3 years, sharp 
variations in image C can be used to plan inten- 
sive monitoring of images A and B to measure the 
response times as a function of wavelength (e.g. 
Kaspi et al. 2007), with the additional advantage 
that the delay between A and B provides redun- 
dancies that protect against weather, the Moon 
and the Sun. The long delays between the im- 
ages also mean that seasonal gaps are completely 
filled, and we can examine the structure function 
of the variability with a densely-sampled, gap-free 
light curve (modulo corrections for microlensing). 
Such data generally do not exist, since most time 
variability data for quasars (other than nearby re- 
verberation mapping targets, e.g. Peterson et al. 
2004) have very sparse sampling (e.g. Hawkins 
2007 on long time scales for a small number of ob- 
jects or Vanden Berk et al. 2004 on shorter time 
scales for many objects). 



In Paper I (jFohlmeister et al. 20071 ) we pre- 
sented three years of optical monitoring data for 
the four brightest images of SDSS J1004+4112 
spanning 1000 days from December 2003 to June 



2006. The fifth quasar image, E, is too faint to 
be detected in our observations. We measured the 
time delay between the A and B image pair to 
be Atba = 38.4 ± 2.0 days. While larger sepa- 
ration lenses tend to have longer time delays, for 
these two images the propagation time difference 
is small, because they form a close image pair (3'.'8) 
from the source lying close to a fold caustic. For 
the more widely separated C and D images we 
could only estimate lower limits on the delays of 
560 and 800 days relative to image B and A. In this 
paper we present the 107 new optical monitoring 
epochs for the 2006/2007 season in §2. When com- 
bined with our previous data we have light curves 
spanning 1250 days that allow us to measure the 
AC delay in §3. In §4 we use the microlensing 
variability of the A/B images to measure the size 
of the quasar accretion disk, and in §5 we measure 
the structure function of the intrinsic variability. 
We discuss the future prospects for exploiting this 
system in §6. 

2. Data 

We monitored SDSS J1004+4112 in the r- 
band during the 2006-2007 season using the Fred 
Lawrence Whipple Observatory (FLWO) 1.2m 
telescope on Mount Hopkins and the MDM 2.4m 
Hiltner Telescope on Kitt Peak. The FLWO ob- 
servations were obtained with Keplercam (0'.'672 
pixels) and the MDM obser vations with RETRO- 
CAM (|Morgan et al. 2005f ) (C/259 pixels). The 
data reduction was carried out as described in 
Paper I. We continued to use the same five stars 
to set the PSF model and the flux scale of each 
epoch and verified that these flux standards con- 
tinue to show no variability. Table 1 presents the 
photometry for the four images in the 2006-2007 
season. 

In Figure[T]we present the resulting light curves 
for images A to D for the period from December 
2003 to June 2007. The average sampling rate 
during the 2006/2007 season is once every third 
day The FLWO data are noisy, so for Figure Q] 
we show a running average of the data (one point 
every five days averaged over ± 7 days) to em- 
phasize the long term trends. Images C and D 
are offset by 0.3 and 1.0 mag, respectively, so that 
they do not overlap with image B in the third and 
fourth season. During this season, images A and 
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Fig. 1. — Light curves of the A, B, C and D images of the quasar SDSS J1004+4112 from December 2003 
to June 2007. Images C and D have been offset by 0.3 and 1.0 mag, respectively, in order to avoid overlap. 
We present a running average of one point every 5 days averaged over ±7 days to emphasize trends and to 
avoid confusion by noise. 



B faded by approximately 0.4 mag with a promi- 
nent feature near the middle of the season, image 
C was relatively constant and image D brightened 
by about 0.4 mag. For the full four seasons, A and 
B have faded by approximately 1 mag, C has re- 
mained relatively constant and D has brightened 
by about 1.5 mag. 

3. The Time Delay 

For the determination of the time delay, we 
use the methods described in Paper I. Our first 
step with the new data was to remeasure the 
A/B delay. The fourth season shows a nice 
feature with maxima in images A and B near 



days 4120 and 4080 respectively, followed by a 
roughly 100 day decline to minima at 4220 (A) 
and 4180 (B) da y s. With the disp ersion method 
(jPelt et al. 1994 iPelt et al. 1991 we measure 
the delay between A and B to be AtsA — 
40.1 ± 3.5 days. For the Kochanek et al. (2006) 
polynomial method we used polynomial orders 
of AT src = 20, 40, 60 and 80 for the source and 
= 1, 2, 3 and 4 for the microlensing variability 
and derived the final estimate using the Bayesian 
weighting of these cases described in Poindexter 
et al. (2007). We found delays of 40.6 ± 1.8, 
40.1 ± 1.8 and 39.8 ± 1.8 (68% confidence regions) 
depending on whether we weighted the changes 
in the number of parameters using the Bayesian 
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information criterion (which strongly penalizes ex- 
tra parameters), the Akaike information criterion 
(which weakly penalizes extra parameters) or no 
penalty for extra parameters. These are consistent 
with our result from Paper I of Atba — 38.4 ± 2.0 
days, but are somewhat more conservative in their 
treatment of the parameterization and the role of 
microlensing. 

In Paper I we derived a lower limit on the BC 
delay of Arcs > 560 days and suggested, based on 
some similarities between the third season for A/B 
with the first season for C, that a delay of order 
700 days was plausible but statistically too weak 
to claim as a measurement. We now see that the 
feature in the second season for image C strongly 
matches the feature we observe in the new season 
for A and B. Using the dispersion spectra method 
(Pelt et al. 1994, 1996), we find At ca = 822 ± 7 
days and Atcb — 780 ± 6 days where the CA de- 
lay is slightly less accurate because the CA over- 
lap is slightly less than the CB overlap due to the 
alignment of the light curves relative to the sea- 
sonal gaps. The three delays are mutually con- 
sistent since Atqb = Atqa — Atba — 782 ± 7 
days. For the polynomial method analysis we si- 
multaneously fit A, B and C holding the A/B delay 
fixed to 40.6 days to find CA delays of 821. 6± 2.1, 
823.0 ± 2.1 and 820.2 ±2.1 days for the three 
weighting methods, respectively. Image D should 
lag the other three images, and we see no feature 
in the light curve of image D that can be matched 
to the first season of images A/B. The lower limit 
on the time delay between images A and D is now 
At da > 1250 days (3.4 years). 

We modeled the lens using the same approach 
as in Paper I, where we include the central cD 
galaxy, and NFW halo for the cluster dark mat- 
ter and 12 pseudo-Jaffe models corresponding to 
cluster galaxies (we added an extra component at 
(x,y) = (31'.'0, 4'.'0) relative to quasar image A in 
an effort to reduce the overall shear). The fits 
were carried out using lensmodel (Keeton 2001) 
and while adequate they are not satisfactory - it 
is very difficult to find solutions with no additional 
quasar images created by the galaxies, and check- 
ing for the extra images makes the procedure ex- 
traordinarily slow. At present we lack the ability 
to model this system in detail (including uncer- 
tainties) at the precision of the constraints, while 
simplified models that ignore the galaxies are inca- 



pable of fitting the data at all. The model predicts 
an AD delay of order 2000 days (5.5 years), which 
is consistent with our current lower bound. 

4. Microlensing and the Size of the Quasar 
Accretion Disk 

The residuals of the A and B light curves 
(see Fig. [2]) clearly indicate that microlensing is 
present. After correcting for the time delay, the 
mean magnitude differences between A and B for 
the four seasons are 0.460 ± 0.005, 0.283 ± 0.007, 
0.339 ± 0.005 and 0.381 ± 0.007 mag. For the 
two seasons overlapping with C we find mean 
magnitude differences, seasonal gradients and sec- 
ond derivatives of 0.590 ± 0.010 mag, -0.04 ± 
0.02 mag/year and 0.29 ± 0.09 mag/year 2 for 
C relative to A and 0.368 ± 0.005 mag, 0.05 ± 
0.01 mag/year and 0.18 ± 0.04 mag/year 2 for B 
relative to A. Fig. [2] shows the superposition of 
the phased A, B and C light curves and the dif- 
ferences between them that are the signature of 
microlensing. 

We modeled the microlensing for images A/B 
using the Bayesian Monte Carlo method of 
Kochanek (2004). We used the microlensing pa- 
rameters of our (adequate) lens model, with con- 
vergence k and shear 7 values of k — 0.48 and 
7 = 0.57 for A and k = 0.47 and 7 = 0.39 for 
B. We allowed the surface density in stars k* 
to vary from 10% to 100% of k increments of 
10%. We used a microlens mass function with 
dn/dM oc M 3 with a dynamic range in mass 
of a factor of 50 that approximates the Galac- 
tic disk mass function of Gould (2000). We 
generated 4096 x 4096 pixel magnification pat- 
terns with an outer scale of 20 (Re) where (Re) 
is the Einstein radius at the mean stellar mass 
(M). We modeled the disk as a face-on, thin 
disk ( Shakura &: Svunvaev 1973h neglecting the 
central temperature depression and relativistic ef- 
fects. We measure the disk size R\ as the point 
where the disk temperature matches the rest- 
frame energy of our monitoring band, kT\ — hc/X, 
where A ~ 2300A for the r-band at the source red- 
shift (see Morgan et al. 2007). The half-light 
radius i?i/ 2 = 2.44i?>, should be used to compare 
to any other disk model, since Mortonson et al. 
(2005) have shown that the half-light radius de- 
pends little on the surface brightness profile of 
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the model. We made four realizations of each of 
the 10 microlensing models and drew 2 x 10 5 trial 
light curves for each of the 40 cases so that we 
would have a reasonable statistical sampling of 
light curves that fit the data well. We found that 

R . = 10 14.8±0.3^^ (1) 

2300A h 70 y cost 

for a disk inclination angle i, whether or not we use 
a prior on the mean microlens mass of O.IMq < 
(M) < M Q . 

From the Mgll emission line width/black hole 
mass calibration of Kollmeier et al. (2006), the 
spectrum of image C from Richards et al. (2004), 
and a magnification-corrected HST /-band mag- 
nitude of 20.9 ±0.4, we estimate a black hole mass 
of log Mbh/M q = 8.4 ± 0.2. Fig. [3] compares the 
disk size estimate to the characteristic scales of 
such a black hole. 

5. The Structure Function 

The quasar structure function can be used as a 
tool to characterize quasar variability independent 
of short-timescale monitoring gaps and to compare 
with theoretical models of quasar variability (e.g. 
Kawaguchi et al. 1998). The structure function 



S(r) = ^^£[™fe)-™(*i)] 2 (2) 

is the variance in the magnitude as a function 
of the time r = tj — U between measurements 
where m(tj) and N(t) is the number of epochs 
at that time lag. For SDSS J1004+4112 we can 
determine the structure function over a moderate 
time range and with a dense sampling rate and no 
seasonal gaps if we use the time-delay corrected 
quasar light curves for images B and C. These 
cover a time-baseline of 2065 days (5.7 years) in 
the observers frame, corresponding to a maximum 
rest-frame time lag at z s = 1.734 of 755 days. 
For the very different behavior of the image D 
light curve, which could not yet be time-delay con- 
nected to the other images, we compute the struc- 
ture function independently for rest frame time 
lags up to 470 days. As in Vanden Berk et al. 
2004 we fit the form of the structure function with 
a power law. The value for the power law index 
7 = 0.35 ± 0.02 for the combined image B and 
C light curves is consistent with that derived for 
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Fig. 2. — The image A, B and C light curves 
in their overlap region after shifting by the time 
delays. The data are binned in one week inter- 
vals. The lower box shows the residual magnitudes 
shifted by the offset between the images, revealing 
microlensing variability of order 0.15 mag. The 
light curve of image B was chosen to have con- 
stant flux because it has the most overlap with 
the over two. 

the SDSS quasar sample. For image D we find 
a similar slope of 7 = 0.39 ± 0.03, as expected 
from the light curve. Time-delay connecting the 
image A, B and C lightcurves by subtracting the 
estimated microlensing variability in the overlap 
region of the lightcurves gives a restframe record 
of the intrinisic quasar variability over 500 days. 
The slope of the structure function for the source 
light curve 7 S = 0.45 ± 0.03 is steeper than for the 
observed non-microlensing corrected curves. 
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Fig. 3. — Probability distribution for the accre- 
tion disk size at 2300A assuming the mean disk 
inclination (cos(i) = 1/2). The solid (dashed) 
curves show the distribution without (with) the 
prior on the average microlens mass. The verti- 
cal lines show the Schwarzschild radius and the 
expected size for thin disk theory assuming the 
black hole mass estimated from the Mgll emis- 
sion line width. The expectation from thin disk 
theory assumes the disk is radiating at Edding- 
ton (L/Le = 1) with efficiency r\ = 0.1 where 
L = r]Mc 2 . 

6. Summary and Conclusions 

We present a fourth season of monitoring data 
for the four bright images of the five image gravita- 
tional lens system SDSS J1004+4112. We confirm 
our previous estimate for the time delay between 
the merging A/B pair, finding that B leads A 
by 40.6 ± 1.8 days. We measure the delay for 
image C for the first time, finding that it leads 
image A by 821.6 ± 2.1 days. We note that 
this is nearly twice the longest previously mea- 
sured delay (the 417 da y delay in Q0957+5 61 
(|Schild k Thomson 19951 : iKundic et al. 19971) ). 
We find a lower bound that D lags A by more 
than approximately 1250 days. Our current mass 
model predicts that D lags A by approximately 
2000 days, which is consistent with the present 
limit. The fractional uncertainties in the AB delay 



are still dominated by sampling and microlensing, 
while the fractional uncertainties in the AC delay 
are dominated by cosmic variance due to density 
fluctuations along the line of sight rather than our 
measurement uncertainties of 0.3% (e.g. Barkana 
1996). 

A detailed model of this system, including the 
constraints from the multiply imaged, higher red- 
shift ar cs (Sharon et al. 2005), the X-ray m easure- 
ments ( Ota et al. 20061 : Lamer et al. 2006 ) and a 
detailed understanding of the uncertainties will 
be a challenge. We lack a completely satisfac- 
tory model for the system at present, in the sense 
that the modeling process is extraordinarily slow 
due to the ability of the gravitational potentials 
associated with the cluster member galaxies to 
generate additional but undetected images of the 
quasar, making it impossible to carry out a reli- 
able model survey. The record of models for this 
system is discouraging. As wc noted in Paper I, all 
three model studies (Oguri et al. 2004; Williams 
& Saha 2004; Kawano & Oguri 2006) generically 
predicted shorter AB delays than the observed 40 
days, and that this could be plausibly explained 
by the absence of substructure (i.e. galaxies) in 
the potential models. The longer AB-C and AB- 
D delays should be less sensitive to substructure. 
Oguri et al. (2004) do not include an estimate of 
the AB-C delays and have A-D delays consistent 
with our present limits. The range of B-C de- 
lays in Williams & Saha (2004) is consistent with 
our measurement of 820 days, but they predict 
AD delays shorter than our current lower bound 
of 1250 days. Kawano & Oguri (2006) predict a 
range for the longer delays over a broad range of 
mass distributions, none of which match our de- 
lays in detail. However, models with sufficiently 
long C-B delays generally have C-D delays long 
enough to agree with our present limits. 

Based on our present mass model we used the 
microlensing between the A and B images to make 
an estimate of the size of the quasar accretion disk 
at 2300A in the quasar rest frame. If we convert 
this to the expected size at 2500A assuming the 
R\ oc A 4 / 3 scaling for a thin disk and assume the 
mean disk inclination cos(i) = 1/2 the scale on 
which the disk temperature matches the photon 
energy is R ° = 10 15 0±0 - 3 cm. Comparisons to 

OJ 2500A 

other disk models should use the half-light radius 
which is 2.44 times larger. Based on the quasar 
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Mgll emission line width we estimate that the 
black hole mass is 1O 8 4±O ' 2 M0. For this mass, the 
microlensing accretion disk size-black hole mass 
correlation found by Morgan et al. (2007) pre- 
dicts that R - = 10 153 cm, which is in broad 

2500A ' 

agreement with the measurement. Further obser- 
vations, the inclusion of additional images, and 
monitoring in multiple bands should improve these 
measurements and potentially allow us to deter- 
mine the mean surface density in stars near the 
images k* and their average mass (M) . Similarly, 
the ability to construct continuous light curves of 
the intrinsic variability and to use image C to pro- 
vide early warning of sharp flux changes that can 
then be intensively monitored in images A and B 
may make this system a good candidate for apply- 
ing reverberation mapping techniques to a mas- 
sive, luminous quasar. At present, we already see 
that the system has a structure function typical of 
quasars. 
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Table 1 

Light Curves for SDSS J1004+4112 * 
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.673 


0, 


,54 


3 


.700±0 


.038 


4. 


016±0 


.050 


4, 


,319±0. 


,067 


4, 


,081±0. 


.054 


FLWO 


Keplercam 


4238 


.688 


0. 


,58 


3 


,713±0 


.042 


3. 


976±0 


.052 


4, 


,485±0. 


,084 


4, 


114±0. 


.061 


FLWO 


Keplercam 


4239 


.745 





,45 


3 


.682±0 


.050 


4. 


147±0 


.075 


4, 


,380±0. 


,094 


4, 


177±0. 


.080 


FLWO 


Keplercam 


4240 


.707 


0. 


,75 


3 


,721±0 


.038 


3. 


985±0 


.048 


4, 


,491±0, 


,076 


4, 


,039±0, 


.052 


FLWO 


Keplercam 


4245 


.660 


0. 


,88 


3 


.715±0 


.053 


4. 


,039±0 


.070 


4, 


,464±0, 


,105 


4, 


,073±0. 


.074 


FLWO 


Keplercam 


4246 


.714 


0, 


.34 


3 


.662±0 


.076 


3. 


887±0 


.092 


4, 


,249±0. 


133 


4, 


,149±0, 


.121 


FLWO 


Keplercam 


4247 


.702 


2. 


.05 


3 


.788±0 


.086 


4. 


,055±0 


.109 


4, 


,247±0. 


,133 


4, 


,130±0. 


,120 


FLWO 


Keplercam 


4248 


.668 


0. 


.77 


3 


.759±0 


.069 


4. 


,235±0 


.104 


4, 


,763±0. 


,169 


4, 


,233±0, 


,107 


FLWO 


Keplercam 


4249 


.684 


1. 


.54 


3 


.850±0 


.081 


4. 


358±0 


.126 


4, 


377±0. 


,133 


4, 


,198±0 


,113 


FLWO 


Keplercam 


4250 


.685 


0. 


.87 


3 


.675±0 


.062 


4. 


109±0 


.091 


4, 


,641±0. 


,147 


4, 


,164±0. 


,098 


FLWO 


Keplercam 


4252 


.675 


1, 


.74 


3 


,531±0 


.067 


4. 


,010±0 


.103 


4, 


,445±0. 


,153 


4, 


.085±0 


,113 


FLWO 


Keplercam 



Table 1 — Continued 



HJD 


X 2 /N dof 




Image . 


A 




Image 


B 




Image C 




Image D 


Observatory 


Detector 


4254. 


.660 





.55 


Q 
O 


.DOZitU 


UoO 


A 

4. 


.U ( ( ±u 


ri/io 
.U4y 


A 
4 


,426±0. 


.069 


A 
4. 


.068±0 


.051 


r JjVV W 


Keplercam 


4255. 


.667 


1 


.34 


o 

O 


(zki _i_n 


,Uo4 


A 
4 


.Uol_tU 


H/l Q 
.U4y 


A 
4 


,431±0 


.068 


A 
4 


.119±0 


.052 


T?T WO 


Keplercam 


4258. 


,675 





.50 


o 
o 


.DZZztU. 


,U4U 




i n/i-i-n 


.uou 


4 


.359±0 


.076 


A 
4 


,144±0 


.064 


r IjVV w 


Kcplcrcam 


4260. 


.697 





.73 


o 
o 


(30Q_|_ri 

.Ooy-tU. 


,UoO 


A 
4 


.U4D_t:U 


.uou 


A 
4 


.395±0 


.069 


A 
4 


.062±0 


.052 


T?T WO 
r IjVV w 


Keplercam 


4261. 


,716 





.51 


o 
o 


.OZO-tU. 


,U4U 


4 


1 1 /i _i_n 


.uou 


4 


.422±0 


.082 


A 
4 


.065±0 


.060 


T?T WO 
r JjVV \J 


Keplercam 


4263. 


.668 





.48 


6 


ftKi _i_n 
.D01±U. 


,uou 


A 


.ZolztU 


.Uoo 


4 


314±0 


.092 


A 
4 


,117±0 


.077 


r JjVVvJ 


Keplercam 


4264. 


,654 





.69 


Q 
o 






-i . 


nqfi-l-n 




A 
-± 


.467±0 


.096 


A 
-± . 


,032±0 


.067 


FT WD 

X ±J V V V t 


XYL.JJ1LX Lcllll 


4265. 


,685 





.48 


3 


,693±0. 


,046 


4 


,141±0. 


.068 


4 


,301±0 


.081 


4 


,114±0. 


.069 


FLWO 


Keplercam 


4266. 


,700 





.73 


3 


,531±0. 


,068 


4 


.243±0 


.129 


4 


,148±0 


.124 


4 


,253±0 


.135 


FLWO 


Keplercam 


4269. 


,687 


1 


.27 


3 


,704±0. 


,083 


4 


,117±0 


.120 


4 


,426±0 


.159 


3 


.972±0 


.111 


FLWO 


Keplercam 


4271. 


,674 





.44 


3 


,742±0 


,077 


3 


.960±0 


.093 


4 


317±0 


.133 


4. 


,037±0 


.104 


FLWO 


Keplercam 


4276. 


,677 





.36 


3 


,720±0. 


,077 


4 


,101±0 


.107 


4 


.360±0 


.140 


4 


,122±0 


.113 


FLWO 


Keplercam 


4277. 


,649 





.25 


3 


,695±0. 


,167 


3 


.986±0 


.213 


4 


.195±0 


.266 


4 


.277±0 


.282 


FLWO 


Keplercam 


4278. 


,668 





.32 


3 


,615±0. 


.079 


4 


.099±0 


.121 


4 


,576±0 


.190 


4 


.002±0 


.116 


FLWO 


Keplercam 



Note. — The Heliocentric Julian Days (HJD) column gives the date of the observation relative to HJD= 
2450000. The \ 2 /Ndof column indicates how well our photometric model fit the imaging data. When x 2 > ^dof 
we rescale the photometric errors presented in this Table by (x 2 /Ndof) 1 ' 2 before carrying out the time delay 
analysis to reduce the weight of images that were fit poorly. The image magnitudes are relative to the comparison 
stars (see text). The magnitudes enclosed in parentheses are not used in the time delay estimates. 



